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Abstract 

In this paper a new approach to solving the Ising-Onsager problem in external 

magnetic field is investigated. The expression for free energy on one Ising spin 

in external field both for the two dimensional and threedimensional Ising model 

with interaction of the nearest neighbors are derived. The representations 

of free energy being expressed by multidimensional integrals of Gauss type 

with the appropriate dimensionality are shown. Possibility of calculating the 

integrals and the critical indices on the base of the derived representations for 

free energy is investigated. 
PACS number(s): 05.50. +q 
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INTRODUCTION 



It is well known that the Ising-Onsager problem JT|, in external magnetic field has 
not been solved till now, despite intensive efforts of a few generations of physicists and 
mathematicians. By the problem we mean the exact calculation of the statistical sum for 
the Ising model with interaction of the nearest neighbors both in the two dimensional case 
(in finite external field) and in the threedimensional case (in external field as well as without 
such a field). Therefore, we do not intend to present the variety of approximate methods 
and approaches to solving the Ising-Onsager problem. Detailed discussion of these matters 
could be found in numerous well known papers and monographs. We note here only the 
paper by C.N. Yang ||, where the problem was investigated for the case of infinitely small 
external field in two dimensions. In our opinion, the efforts to find an exact solution to 
the Ising-Onsager problem is of great interest till now. The reason is, the Ising models are 
connected with numerous other models of statistical mechanics and quantum field theory 
@] — J7|. Therefore the present paper we treat as the next step in the direction described 
above. 



I. THE PARTITION FUNCTION 

Here we consider firstly the two dimensional case (d = 2, H ^ 0), and then we will report 
on the results for the threedimensional case (d — 3, H ^ 0) without going into details of 
derivations. Let us consider a rectangular lattice consisting of M columns and N lines, in 
nodes of which are given variables a nm which take values ±1. 

These variables will be called below "spins". The collective index nm numbers nodes of 
the lattice; n numbers a line, m numbers a column. The Ising model with nearest neighbors 
interaction is given by the following form of the Hamiltonian: 
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which takes into account possible anisotropy of the interaction between nearest neighbors and 
also interaction of spins a nm with external field H, which is directed "upwards" (cr nm = +1). 
The investigated problem consists of calculation of the statistical sum for the system: 

Z(h)= £ ... £ exp(-/3W) 

<TU=±1 "\2VM=±1 



£ exp 

(o"nm=±l) 



NM 



n,m=l 



[1.2) 



where 



#i,2 = PJifi, h = PH, (3 = l/k B T. 



;i.3) 



Typically periodic boundary conditions on variables o nm are imposed and we will assume 
this everywhere below. Let us note that the statistical sum (1.2) is symmetric with respect 
to the change h — > —h. 

As is known ||, the statistical sum (1.2) can be represented in the form of the trace of 
the T-operator (T-transfer matrix): 



Z(h) = Tr(T) M = Tr \(2 smh2K 1 ) {N/2) T 1 T 2 T h 



M 



;i-4) 



where the matrices T\ 2 h are of the form: 



^i = exp(^fx) 



;i.5) 



N 



T 2 = exp K 2 £ « +1 , a z N+1 = a*, 



n=l 



T h = exp £ a^j 



and Kl and K\ are connected by the relations: 



;i.6) 



;i.t) 



exp(-2Ki) = tanh(ir*)> sinh(2iTi) sinh(2ir*) = 1. 



;i.s) 
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In the formulae (1.5) — (1.7) the quantities (<j'n ,y ' z ), (n = 1,2, ...N) - are well known from 
quantum mechanics 2 JV -dimensional matrices: 



fjX.V.z =1^1 



T x,y,z 



where a x,y ' z — twodimensional spin Pauli matrices: 



v 1 V 



i 



1, (JV factors), 



<1 0^ 



a 



v° - 1 / 



(1.9) 



satisfying the standard transposition relations: 



a k f = 1, o-V + crV* = 0, (k,j) = x,y,z; ;a x a y = ia z , 



;i.io) 



For example, spin matrices for the n-th electron in the system consisting of N nonrelativisitc 
electrons are exactly the matrices a x ' y,z . It is known that for n ^ n spin matrices a x,y ' z 
commute and for any given particular n they satisfy formally relations (1.10). It follows 
from this that matrices T 2 and 7\, (1.6 — 7) commute but they do not commute with the 
matrix Ti, (1.5), i.e. 



(1.11) 



It follows from (1.11) that under Tr(...) we can write for the statistical sum (1.4) the 
expression: 



Z{h) = Tr{jTiTl /2 T 2 Tl ,2 ) M = Tr('yT^ 2 T 1 T^ 2 T 2 ) M = Tr(P) M , 



(1.12) 



P = {2sinh2K 1 ) N/2 Tl /2 T 1 Tl /2 T 2 , 



1.13) 



where we used the identity Tr(AB) = Tr(BA). 



Now we will consider in more details the matrix U = T^T^ 2 ; since the matrices a: 



and a z ,, commute for n ^ n', we can write the matrix U in the form: 

N N 

U = T]! 2 T{T]! 2 = I] e W2K e *r«* e (V2K = J] [/ n . 



;i.i4) 



n=l 



n=l 



Further the matrix U n can be represented in the following form: 

U n = e (/l/2),J ™(coshK* + <sinhK*)e ( ' l/2)CT " 



exp 



/ cosh K{ sinh h sinh K\ x 

U ~i °V H ~i a n 

\ sinhw smhw , 



1.15) 



where ui - is a positive root of the equation: 



cosh uj = cosh X* cosh h, 



In determination of the expression (1.15) we used the identity: 



exp (/it) = cosh fx + t sinh /i, t 2 = 1. 



1.16) 



1.17) 



It is easy to see that for ft = we obtain from (1.12) the standard expression for the 
statistical sum Z for the twodimensional Ising model without external field 0. [J . 

A. The ID Ising model 



Relatively easy one can show that (1.12) becomes: 



Z(h) = Tr 



N 



(2 sinh 2K X ) N ' 2 J] U n T 2 



n=l 



M 



1.18) 



where U n is given by the formula (1.15), and it describes correctly the transition to the 
onedimensional Ising model both in the constant Ki, and in the constant K 2 - Indeed, if we 
take K 2 = and N = 1, and this procedure corresponds to neglecting summation over n we 
obtain: 



Z x {h) = Tr [(2sinh2ir 1 ) 1/2 f/ c 



M 



1.19) 



where the matrix Uo is given by (1.15), where the index n was omitted. Eigenvalues of the 
matrix Uq can be easily obtained: 



A = exp(±u;) 
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where u — is the positive root for the equation (1.16). As a result we obtain the following 
formula describing free energy on one spin in the thermodynamic limit: 



f(h) = -- lim — lnZiO) = -ilnfe^cosh/i + fe^sinh^ + e-^ 1 ) 1 / 2 



;i.2o) 



i.e. the known classical expression ||T|. 

Transition to the onedimensional Ising model limit in the constant K\ could be done by 
taking K\ — > 0, and M — 1, i.e. we neglect summation over the index m and go to the limit 
K\ — > 0. As a result we get from (1.18): 

JV 



Z x {h) = Tr 



n(i +ot 2 t, 



.n=l 



:i.2i^ 



where we used the transition to the limit: 



Kmj2 sinh 2K 1 ) 1/2 &cp{ K K) 



where we took into account the relation (1.8) between K\ and K*. It is reasonable to stress 
here that the factors (1 + a*), entering the expression (1.21), are simply necessary to get 
the correct result. For the aim of calculating the trace (1.21), it is convenient to go to the 
fermion representation [§],§[. Omitting some calculations we write for Zi(h), (1.21) the 
expression: 



Z 1 (h) = Tr(DT±T h ) 



[1.22) 



where the operators D, and T^, expressed in terms of Fermi creation and annihilation 
operators (c+, c n ) are of the form: 

N 



;i.23) 



n=l 



it = exp 



N 



K 2 zZ( C t ~ C n)(4+1 + C n +l] 



n=l 



;i.24) 



( N 


n— 1 




I h ^2 exp 


Z7T C+C p 




[ n=l 


p=l 





;i.25) 
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In the formula (1-24) the sign (+) is related to states that are even with respect to the 
operator of the complete number of particles (N = J2n=i c t c n)> to which there correspond 
anticyclic boundary conditions, and the sign (— ) to the odd states, to which there correspond 
cyclic boundary conditions. It is easy to see that because of the multiplicative character of 
the operator D, (1.23), all diagonal matrix elements in (1.22) vanish with the exception of 
the vacuum-vacuum matrix element, i.e.: 



Z x (h) = 2 N < | (TfT h ) | >, 



1.26) 



where the operators T 2 and Th are defined by (1.24 — 25). Then, "acting" with the operator 
Th on the vacuum state | >, and using the Hausdorff-Baker formula ||10|| , (a, (3 = const): 

exp(aa;) exp(Py) = exp(ax + fly + (a/3/2) [x, y]~), 



[x,z]- = [y,z]- = 0, z = [x,y]-, 

the operator Th, (1.25) can be reduced to the "effective" form (in the sense of action on 
|0>): 



T h = cosh (h) exp 



N N-n 

tanh 2 hJ^Yl c n c n 

n=l p=l 



+P 



;i.27) 



When developing the expression (1.27) we have taken into accaunt the fact, that the diagonal 
matrix elements of the odd number of Fermi operators are equal to zero. 
Finally, going to the momentum representation: 

- - exp( -" /4, E^, 



and computing the matrix element for a fixed q after some not complicated transformations 
we arrive in the case of even states at the expression for the statistical sum [Z^(/i)] (1.26): 

Z+(h) = [2cosh(/i)] JV J] [cosh 2fsT 2 - sinh 2fsT 2 cos q + a 2 sinh2AT 2 (l + cos q)) 

0<q<TT 

N r O -,1/2 



[2cosh(/i)cosh J fs: 2 ] A ' JJ 



n=l 



1 + z 2 + 2z 2 z — 2z 2 (l — z) cos 



27m, 



1.28) 



where z 2 = tanhi^ and z = a 2 = tanh 2 h. In the case of odd states it is easy to show that 
the sum Z-y {h), is equal to Z-y [K] — IZ^iK). Finally, we obtain in the thermodynamic limit 
again the formula (1.20), (M — > N, Ky —> K 2 ) for free energy on one spin. 

The ID Ising model was discassed here with so many details because it was unexpectedly 
found, that Zy(h) such as represented in (1.28) can be applied in graph theory. Namely, 
using the representation (1.28) one can calculate the generating function for the Hamilton 
cycles on the simple square lattice (N x M) , || . 



B. The 2D Ising model 



For further purposes it is convenient to represent the matrix UT 2 , entering the formula 
(1.18) and where U is defined by (1.14) in the form of a simple product of matrices P n such 
that their diagonalization is relatively easy. The reason is that, as is known from — |§, 
to calculate free energy on one spin in the thermodynamic limit it is sufficient to find the 
maximal eigenvalue of the matrix UT 2 , which is 2 N x 2 N dimensional. First of all we note 
that the matrix U (1-14) can be represented in the form of a simple product of matrices Uq. 



N 



U = II U n = u o ® U <8> ... <8> U , N - factors, 



;i.29) 



71=1 



where the matrix Uq is defined by the formula (1.15), in which one should skip the index 
n. In order to represent the matrix T 2 (1.6) in the form of a simple product we will use the 
well known identity |nj 11 j: 



1 f°° 

exp(A 2 ) = -— / ex V (-e + 2A£)d£, 

7T ' J-oo 



;i.3o) 



where A — is a bounded operator (matrix). Writing exp(K 2 a^a^ +1 ) in the form 

K 2 



exp(fs:2« + i) = exp 



we can represent the matrix T 2 in the form: 



71 



-NK 2 
N/2 



exp 



A? 



~2 + a n+l, 



N 



K 2 



E^ + (2^ 2 ) 1/2 E(^ + e, + i)< +1 



71=1 



n=l 



A' 



;i.3i) 



n=l 



where a z N+1 = a\ and £at+i = £1. After writing the matrix T 2 , (1.31) this way we can 
represent it in the form of a simple product of matrices exp[(2K 2 ) 1 ^ 2 (^ n + £n+i)o~ z ] inside the 
integral: 

N N 



I] exp[(2Jf 2 ) 1/2 (e„ + £n + iK +1 ] = J] ® exp[(2Jf 2 ) 1/2 (£„ + £ n+l)(7 *], 



;i.32) 



n=l n=l 

where on the right hand side of the formula there is a simple product of 2 x 2 matrices. 
Next, we can write the matrix UT 2) using (1.29) and (1.31 — 32), in the form: 



e -NK 2 f00 roo N r N i r N 

UT * = / •- / II ex P " £ ^ II ® exp[(2K 2 ) 1 /2 (a + 

7T / J-oo J-oo„ =1 L n=l J Ln=l 



, (1.33) 



where we included the constant matrix U under the integral and we used the known theorem 
on simple product of matrices: 



(Ai ® A 2 <g> ...)(Bi <g> B 2 



(AiBi) ® (A 2 B 2 ) 



Using the expression (1.33) enables calculation of all 2 N eigenvalues of the matrix UT 2 . 
Eigenvalues of the matrix U exp[a(£ n + £„+i)o" 2 ] can be easily calculated and are equal to: 



A ± (n,n + 1) = e ± ^"' n+1 \ 



where cu(n, n + 1) is defined as a positive root of the equation: 



cosh[u;(n, n + 1)] = cosh(^) cosh[/i + a(£„ + £„+i)], « = (2i^ 2 ) 



1.34) 



1/2 



;i.35) 



In the diagonal representation the matrix V under the integral (1.33) can be represented in 
the form: 

/\ + , ^ « \ r 



V = 



A 7 " 



[J ®S(n,n + l) 

.n=l 



N 

n 



J n=l 



IJ ®S"(ra>" + 1 ) 

Ln=l 



, (1.36) 



\ + (n,n + l) 
^ \-(n,n + l) J 

where S(n,n + l)S'(n,n + 1) = 1, and X ± (n, n + 1) are defined above by (1.34). From this 
it follows that the eigenvalues Aj of the matrix V are equal to: 



A, = A ± (l, 2)A ± (3, 4)...A ± (A, 1), (j = 1, 2, 3, 2 JV ), 



(1.37) 



where to each j there corresponds a combination of (+) and (— ) eigenvalues A ± (n, n + 1). 
Finally we can express the statistical sum (1.18) by the formula: 



Z(h) = Tr 



-NK 2 ,00 ,00 / N - W 



poo poo " I \ 

(2smh2^ /2 ^ / ..../ n*n«p -E^r 



;i.38) 



where the matrix V is given by (1.36). 



II. THE FREE ENERGY 



As we mentioned above, the free energy on one spin in the thermodynamic limit can be 
expressed by the maximal eigenvalue of the matrix UT 2 , entering (1.18). We menaged to 
express this matrix in the form of an A/"- type integral (1.38), where the matrix V is defined 
by (1.36), and all the matrix elements of V are positive. On the other hand, in accordance 
with the known Frobenius- Perron theorem, the matrix B, with all matrix elements positive 
has its maximal eigenvalue nondegenerate. Let us assign to the maximal eigenvalue of the 
matrix V a letter A max . In accordance with our definition of the eigenvalues X ± (n, n + 1), 
using (1.37) we obtain the following expression for A max : 

N N r N 

A max = n A + (n, n + 1) = I] ^ l ' n+l) = exp £ cu(n, n+1) 

n=l n=l Ln=l 

where u(n,n+ 1) is defined as a positive root of the equation (1.35), and 



(2.1) 



A- m ax > Aj, (j = 1,2,...,2 W ). 

Further we denote the eigenvalues of the matrix UT2 by Aj. Taking into account the dimen- 
sion of the matrix UT 2 , which is equal to 2^, we can write on the base of the relation (1.18) 
obvious inequalities: 

< Z(h) < 2 N A^ ax , (2.2) 

where A max - - is the maximal eigenvalue in the set Aj, to which we included also the 
constant factor (2 sinh2A'i) Ar / 2 . Taking the logarithm (2.2) of this expression and dividing 
by the nodes number NM, we arrive at the next system of inequalities: 
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^ ln(A maa; ) < In Z(h) < 1 ln(A_) + ^ In 2, (2.3) 

in which the expression in the middle representations free energy on the node if we neglect 
the factor — /5 1-1 , where (3 = T is temperature. Going to the limit (N,M) — > oo, we 
obtain the desired formula describing free energy on one spin in the thermodynamic limit: 

where A maa; is, in accordance with (1.38) and (2.1) equal to: 



e ~NK 2 ^oo poo " 

A max = (2sinh2K 1 ) Ar / 2 ^— / .... / J] df n exp 

7T ' J— oo J— oo 

Finally, using the Onsager identity: 

1 

| x |= — / dq\n[2 cosh(x) — 2 cos(g)], (2.6) 
7T Jo 

we obtain the following expression for the function u(n,n + 1) (1.35) 

u(n,n + l) = - r dq\n[2coshK*cosh(h+ {2K 2 ) 1/2 {£ n + £ n+1 )) -2cos(g)], (2.7) 

IT JO 

Expressions (2.4) and (2.5) should describe properly at least the transition to the onedi- 
mensional Ising model. It is easy to show that the transition to the limit K 2 = 0, gives the 
correct result (1.20) for the onedimensional Ising model. The analogous limit taken with 
respect to the constant K\ seems a little bit more complicated and has the form: 



(2.5) 



lim KmaxiKi) = / .... / 11 d£ n exp 

Ki— »U 7T ' J— oo J — oo ^ 



N 



-NK 2 r oo roo N 



n=l 



N 



J\2 cosh[/i + a(£ a + £ n +i)], 



n=l 



(2.8) 



where a is defined above by (1.35). This is an integral of the Gauss type and it could 
be relatively easily calculated. For this purpose we apply the following formal procedure. 
Namely, let us write the expression 2cosh(...), entering the integral (2.8), in the form: 

2cosh[/i + a(£ n + £ n+1 )] = exp[^ n h + afi n {^ n + ^ n+1 )}, (n = 1, 2, N), (2.9) 
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where we introduced a new variable // n of the Ising type. Therefore, we can represent the 
right hand side of the equality (2.8) in the form: 



E 

On=±l) 



-NK 2 



7T 



N/2 



N 



Yl di n exp 



1 n=l 



A? 



^ y £71 



71=1 



A? 



71+1/ 



71=1 



E ex p 

(M»=±l) 



A 



.71=1 



(2.10) 



where we took an integral over the variables £ n , and we imposed on variables fi n cyclic 
boundary conditions (/ijv+i = A*i)- Calculation by standard methods |||7] of the sum 
(2.10), and following substitution of the expression (2.4), gives well known result (1.20). 

Consideration of the expressions (2.4) and (2.5) for free energy of the two dimensional 
Ising model in external field we present in the end of this paper but now we go to the 
threedimensional case. 



III. THE THREE-DIMENSIONAL ISING MODEL 

Hamiltonian for the threedimensional Ising model in external field with nearest neighbors 
interaction we write in the form: 

NMK 

T~l = — ^ {J\0'nmk a 'n,m+\,k + J2 cr nmk cr n+l,mk + ^3 cr nmfc cr nm,fe+l + Ha nm k) , (3-1) 
(n,m,fc)=l 

where the collective index (nmk) numbers nodes of the simple cubic lattice and H is the 
external field. Constants Ij take into account anisotropy of interaction of Ising spins. We 
impose on the variables a n mk, as it is commonly done, periodic boundary conditions. Quan- 
tities N, M and K are node numbers in corresponding directions of a cubic lattice. As is 
known ||, the statistical sum for the threedimensional Ising model can be represented in 
the form of a trace of the A"-th power of the fiber- fiber transfer matrix (R): 

W(h) = Tr(R) K = Tr(T 3 T 2 T 1 T h ) K , (3.2) 

where the matrices T,, (i = 1, 2, 3, h) of dimensions 2 NM x 2 NM are of the form: 
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Ti = exp K x V < m < m+ i , T 2 = exp [ K 2 V < m <7, 



/ ; nm u n,m+l 
nm / 



nm / 



T 3 = (2 sinh 2K 3 ) JVM/2 exp (V* ]T < m ) , T fc = exp (fc £ 

\ nm / \ nm 



0", 



(3.3) 



(3.4) 



Here K t = (3J h (i = 1,2,3); (3 = (l/k B T), T - temperature, h = /3H, and K 3 and are 

connected by relations of type (1.8). In the formulaes (3.3 4) the matrices are Pauli 

matrices, which are defined analogously to (1.9), and have dimensions 2 NM x 2 NM . 

Continuing considerations analogous to these in the twodimensional case we obtain the 
following formula describing free energy on one spin in the thermodynamic limit: 



fi(h) — lim lnA moT , 

^ V ' (3 N,M^oc NM 

where the maximal eigenvalue A max of the matrix R, (3.2) is defined by: 



(3.5) 



(2smh2K 3 ) NM / 2 - 



-M(Jf!+)f 2 ) 



7T 



NM 



OO J —CO 



II dr lnmdinm exp 



NM 



E (-V 2 nm-£ m + ^( n , m )) 



n,m=l 



(3.6) 



and u(n,m) is defined aa the positive root of the equation: 



C0shu(n,m) = COsllfQ C0Sh[/l + 0!i(?7 n+ i, m + ?fo+l,m+l) + Mfn.m+l + £n+l,m+l)], (3.7) 

where a h2 = (2K 1:2 ) 1/2 . 

We impose on integration variables r\ nm and £„ m cyclic boundary conditions, in accor- 
dance with periodic boundary conditions for the former variables: 



a 



N+l,m 



<, m , <m+i = +<i, n(m) = 1, 2, 3, ....N(M). 



(3.8) 



Similarly as in the twodimensional case, the function u(n,m) can be expressed explicitly in 
terms of variables i] nm and £ nm , using for this aim the integral representation given by the 
Onsager identity (2.6). 
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It is a little bit more complicated matter to take a limit with respect to the constant of 
interaction K 3 , although the derived formula is much simpler than the formulae (2.4 — 5), 
Imposition of the limit (K 3 — > 0) in the formulae (3.5 — 7), gives, after some simple transfor- 
mations, the following representation of free energy on one spin of the twodimensional Ising 
model /2(h): 

f2(h) — — — ( lim ) In A max = 

1 1 { e ~ NM ( K l+ K 2) poo roo 

~~R(fj ^ ATM 111 1 ^NM / • • • / IT dVnmdtnm exp [- ]T (^ m + f £ )] 

fj (N,M)^oo MM [ 7T JVM J -00 J -00 nm 

X Yl 2 COSh[/l + ai(?7 n+ i im + 1] n+hm+1 ) + a 2 {£ n ,m+l + fn+l,m+l)] > , (3.9) 
nm ) 

where we used relations of the type (1.8), and a 1}2 are defined above (3.7). The integrals in 
(3.9) are integrals of the Gauss type and, as it is easy to show applying the described above 
formal way of introducing a variable of Ising type \i nm = ±1, can be represented in the form 
(1.2). 

Analogously, one can show rigorously that free energy on one spin for the threedimen- 
sional Ising model can be represented in the form of a multiple integral of the Gauss type: 

I l r e -NMK{K 1 +K 2 +K z ) -oo 
h{h) = -p^Jffi^NMK^ 7-00-7-00 I! d ^mAm k dCnm k 



nmk 



X exp[ ^ (Vnmk + inmk + (nmk) 
n,m,k 



x Yl 2cosh[/i + ai(r] nmk + r] n+limk ) + a 2 (£ nm k + £n, m +i,k) + a 3 (( nmk + ( nm ,k+i)] } , (3.10) 

nmk ) 

where on = {2K^I 2 \ (i = 1,2,3). The formulae (3.10) can be in obvious way generalized to 
describe d-dimensional Ising models but we will not stop on these matters here. 

We will make here a few remarks. First of all, as far as it is known to the author, 
the representations for free energy on one spin in the forms (2.4 — 5) and (3.5 — 10) for 

14 



the two dimensional and threedimensional Ising models, respectively, did not appear in the 
literature. We belive that the known representations (see, e.g. JT2HT1]) are more complicated 



than the ones derived by us. For example, in the paper jnj] was given the following integral 
representation for the Ising model: 

/OO ] 
exp(--bs 2 m - «4J exp(K s n s n +i)ds m , (3.11) 

-°° 1 n,i 

where (m, n, i) — - are vector indices ( we use notation from the paper [Q)- I n the limit 
u — > oo and b — > — oo, (b = —4m) the Ising model is recovered (for this aim one should use 
additionally the fact that to every particular spin there is a factor (u/n) 1 ^ 2 exp(— u)). It 
is clear that the presence of the term in exp(...) of the expression (3.11) considerably 
complicates analysis of this expression. Further, although the formulae (3.9) and (3.10) are 
in a sense obvious, the formulae (2.4 — 5) and (3.5 — 7) are not obvious. The integrals 
(2.5) and (3.6) can be represented as integrals of the "quasi Gauss" type, because the 
functions u(n,n + 1) and u(n,m), described by the relations (1.35) and (3.7), respectively, 
in accordance with the Onsager identity (2.6) are almost "linear" in their arguments (£ n ) and 
{Vnm,£,nm)- This justifies our hopes we could learn now to calculate rigorously the integrals 
(2.5) and (3.6) in case h = using an Ising type variable (/i = ±1), described above. On 
the other hand, it seems to us that for the case (h ^ 0) it is much simpler to deal with the 
expression (2.5), than with the expression (3.9), although it could sound paradoxically. For 
the threedimensional Ising model in external field (h ^ 0) the situation is no longer so clear, 
for the infinitely small field (h ~ 0), similarly as in the two dimensional case, it is easier to 
analyze the expression (3.6), than (3.10), that we shall be show in a next publication. 

IV. CONCLUSIONS 

The derived expressions (2.4 — 5) and (3.5) — (3.10) for free energy on one spin for the 
Ising model can be of some interest, we hope. The reason is we actually should learn how 
to calculate logarithmic asymptotics of multiple integrals of Gauss type, as it could be seen 
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from (2.4) and (3.5). It is known, that there exist a well developed formalism of calculation of 
logarithmic asymptotes for integrals of the Laplace type for the onedimensional as well as for 
the multidimensional cases In the case under consideration the situation is a little bit 

more complicated, because for the integrals of the kinds (2.4) and (3.5) it is not possible to 
transform them to a form of a multiple integral of the Laplace type at least in the framework 
of their classical definition. With the increase of the large parameter (A — > oo) there changes 
also the number of variables, over which one integrates and this needs reformulation of 
the corresponding methods of asymptotic estimation of the considered integrals In 
future publications we intended to investigate in more details the expressions (2.4) and (3.5), 
obtained in this paper and to calculate critical indices for the Ising model. 
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